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In birefringent silica fibers, stable propagation of a monochromatic wave can be inhibited 
by a nonlinear process called vector modulation instability (V-MI) in both dispersion regimes 
(normal or anomalous) This contrasts with the scalar modulation instability (S-MI) 
that does not require birefringence but can only arise in the anomalous dispersion regime 
(at least when second order dispersion dominates^). Two limits, those of weak and strong 
birefringence, are amenable to relatively simple analytical study. 6 These predictions have 
been confirmed experimentally in a number of cases,— particularly in the normal dispersion 
regime. 

The only experimental investigation of V-MI in the anomalous dispersion regime that we 
are aware of is a recent unsuccessful attempt using photonic crystal fibers.— Here we report 
what is to our knowledge the first experimental study of V-MI in strongly birefringent silica 
fibers in the anomalous dispersion regime. We also carry out a very precise comparison 
between experimental results and the predictions of numerical simulations. 

Modulation instabilities (MI) can be induced by classical noise present initially together 
with the pump beam. But MI can also arise spontaneously through amplification of vacuum 
fluctuations. 11 In practice classical input noise and vacuum fluctuations compete for inducing 
MI. The experiment reported here is carried out in the regime where the quantum noise is 
dominant. 

Elsewhere,— we present an unified approach to the problem of scalar and vector MI based 
on the stochastic nonlinear Schrddinger equations (SNLSE) which generalizes the work of 
Ref. This approach is particularly well suited for numerical simulations in complex 
situations where classical noise and vacuum fluctuations act together, where the pump is 
depleted, or where the higher order harmonics of MI appear. In previous work on modulation 
instability, comparison between theory and experiment has generally been limited to noting 
that the frequency at which the maximum gain of the MI occurs is correctly predicted. 
Here we show that there is excellent agreement between the experimental results and the 
numerical integration of the SNLSE. In particular the detailed shape of the output spectrum 
can be predicted in detail, even in the case where higher order harmonics appear (which 
cannot be predicted by linear perturbation theory). To our knowledge this is the first time 
experimental and theoretical studies of MI are compared in such detail. A related work 
is the comparison between theory and experiment for RF noise measurements reported in 
Ref.Q- 
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The experimental setup is reported in Fig. [T] It consists of a Q switched laser (Cobolt 
Tango) that produces pulses at 1536 nm, with a 3.55 ns full-width-at-half-maximum 
(FWHM) duration r and a 2.5 kHz repetition rate /. The Gaussian spectral shape of the 
laser has been characterized using a Fabry-Perot interferometer. The measured 0.214 GHz 
FWHM spectral width is slightly larger than expected in the Fourier transform limit. The 
pump power is adjusted using variable neutral density filters (ND). We measured the in- 
jected mean power P m at the end of the fiber. The peak power Pq is relied to the mean 
power according to the relation: 



A polarizing beam splitter (PBS1) ensures the pump pulse is linearly polarized. A half- wave 
plate is used to tune the angle 9 between the pump polarization direction and the principal 
axes of the fiber. A polarizing beam splitter (PBS2) can be used in order to observe the field 
components polarized along the fast or slow axes separately. Lastly, spectra are recorded 
using an optical spectral analyser (OSA). In our experiment we use the Fibercore HB1250P 
optical fiber. The fiber length L = 51 m, the group-velocity dispersion (GVD) parameter 
fa = —15 ps 2 km^ 1 and the group-velocity mismatch parameter A/3i = 286.1 fs m -1 have 
been measured by independent methods (only significant digits have been indicated). Note 
that the accuracy on the value of (3 2 is poor compared to the standards. This is because 
the interferometric method 15 that we used turned out to be difficult to implement with a 
birefringent fiber. The group-velocity mismatch parameter A/?i is deduced from the walk-off 
between pulses propagating on the principal axes of the fiber. The fiber length L is deduced 
from a measurement of the pulse time of flight. The other important parameters of the fiber, 
and a more accurate estimation of fa, can be inferred from MI spectral peaks positions, as 
explained further. 

Fig. El shows a typical spectrum at the fiber output when the angle 9 is set to 45 degrees. 
The fast and slow polarization components have been separated using PBS2 and their spectra 
recorded successively. The plot clearly exhibits two V-MI peaks at 1511.4 nm and 1561.4 nm 
that are polarized along the fast and slow axes respectively. It also shows S-MI peaks at 
1530.0 nm and 1541.9 nm, with first harmonics. In contrast with V-MI, S-MI is equally 
generated on both principal axes. By polarizing the input field along the fast or slow axes, 
we have observed that V-MI disappears and that the amplitude of the S-MI peaks increases 
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dramatically (figure not shown). 

According to linear perturbation analysis, the angular frequency shifts from the pump of 
the MI peaks are given by 



for S-MI and V-MI peaks respectively. Here, 7 stands for the Kerr nonlinearity parameter of 
the fiber. Fig. El shows the evolution of the spectrum of light emerging from the fiber when 
the pump power is increased. Using Eqs. (J2J) and (j3J), the ratios f^r = 18.740 (rad) THz 
and nrr = 0.2135 (rad) THz 2 W _1 where deduced from these measurement. The first ratio 



and the measured value of A/?i permits to infer that fa = —15.27 ps 2 km , which is 
compatible with the independently measured value. From the second ratio, we deduce that 
7 = 3.26 W" 1 km' 1 . 

The exponential growth of the MI peaks and harmonics is clearly apparent on Fig. |3J 
From these measurements we deduce that the ratio between the maximum gain of the V-MI 
and of the S-MI is 0.67±0.05, in good agreement with the theoretical value 2/3. We also find 
that the ratio between the maximum gain of the 1st harmonic and of the S-MI is 1.88±0.15, 
in good agreement with the theoretical valu e 16117 of 2. 

We now focus on the quantitative comparison between experimental spectral amplitudes 
and those predicted by the SNLSE model for spontaneous (or vacuum-fluctuations induced) 
modulation instabilities. This comparison makes sense because the exact shape of the spec- 
trum, and in particular the relative intensities of the modulation instability peaks and har- 
monics, is very strongly dependent on the initial noise and pump peak power. Experimental 
and computed spectra are plotted together in Fig. HJ In the simulations we used the pa- 
rameters deduced from experimental MI peaks positions (see above), but in order to obtain 
a good fit we had to increase the peak pump power by 5% with respect to that deduced 
from the measurements using Eq. We are not sure of the origin of this discrepancy. 
It could either be due to a systematic error in the measured values of P m , to an error in 
the experimental measure of A fa, to the fact that the experimental pulses are not exactly 
Fourier-transform limited, or to some classical noise photons (for instance due to Raman 
scattering in the fiber) that are added to vacuum fluctuations and slightly speed up the 
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instability. In any case the discrepancy is small enough to confidently conclude that in our 
experiment the MI is mainly induced by vacuum-fluctuations. Indeed with this small ad- 
justment the experimental MI spectra are very well reproduced by numerical integration of 
the SNLSE model. 

In summary, we report what is to our knowledge the first experimental observation of 
spontaneous vector modulation instability in a highly birefringent silica fiber in the anoma- 
lous dispersion regime. The pump power dependence of the detuning of both scalar and 
vector side-bands, as well as their polarizations, agree with linear perturbation theory when 
the pump depletion is small. We have also obtained very good agreement between the exper- 
imental spectra and those obtained by numerical integration of the SNLSE derived from the 
quantum theory. This is to our knowledge the first time that theoretical and experimental 
spectra are compared in such quantitative detail. This very good agreement between the two 
approaches proves that the modulation instability that we observed was truly spontaneous, 
in the sense that it mainly results from the amplification of vacuum-fluctuations. 
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List of Figure Captions 

Fig. 1. Experimental setup. ISO: isolator, ND: variable neutral density filter, PBS: 
polarizing beam splitters, PMF: polarization maintaining fiber, and OSA: optical spectrum 
analyser. 

Fig. 2. Spectra of the fast (black curve) and slow (gray curve) polarization components 
of the light emerging from the fiber. The light is injected with a polarization angle 9 = 45° 
relative to the principal axes of the fiber. The mean power is approximatively 1 mW. The 
resolution bandwidth is 1 nm. 

Fig. 3. Spectrum of the output field for increasing pump mean power, respectively 
1.065 mW, 1.105 mW, 1.145 mW, 1.165 mW, and 1.265 mW. The light is injected with 
a polarization angle 9 = 45° relative to the principal axes of the fiber. The resolution 
bandwidth is 0.1 nm. 

Fig. 4. Comparison between experimental spectra (black curves) and numerical inte- 
gration of the SNLSE (grey curves). The numerical results are noisy because only one 
realization of the stochastic method has been computed for each curve. The flat parts on 
the experimental spectra correspond to the sensibility limit of the OSA. In the simulations 
the pump pulse is assumed to be Fourier transform limited. The simulation parameters are 
A = 1536 nm, r = 3.55 ns, L = 51 m, 7 = 3.26 W^km" 1 , (3 2 = -15.26 ps^m" 1 and 
A(3i = 286.1 fs m -1 . The peak powers P corresponding to simulations (experiments) are 
(a) 132.5 W (125.6 W), (b) 127.2 W (121.4 W), and (c) 116.6 W (112.9 W). The numerical 
results have been convolved with a response function to take into account the resolution of 
the OSA. 
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Fig. 1. Experimental setup. ISO: isolator, ND: variable neutral density filter, PBS: po- 
larizing beam splitters, PMF: polarization maintaining fiber, and OSA: optical spectrum 
analyser. Amans57847Fl.eps 
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Fig. 2. Spectra of the fast (black curve) and slow (gray curve) polarization components of 
the light emerging from the fiber. The light is injected with a polarization angle 9 = 45° 
relative to the principal axes of the fiber. The mean power is approximatively 1 mW. The 
resolution bandwidth is 1 nm. Amans57847F2.eps 
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Fig. 3. Spectrum of the output field for increasing pump mean power, respectively 1.065 mW, 
1.105 mW, 1.145 mW, 1.165 mW, and 1.265 mW. The light is injected with a polarization 
angle 9 = 45° relative to the principal axes of the fiber. The resolution bandwidth is 0.1 nm. 
Amans5 7847F3 . eps 
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Fig. 4. Comparison between experimental spectra (black curves) and numerical integration 
of the SNLSE (grey curves). The numerical results are noisy because only one realization of 
the stochastic method has been computed for each curve. The flat parts on the experimental 
spectra correspond to the sensibility limit of the OSA. In the simulations the pump pulse 
is assumed to be Fourier transform limited. The simulation parameters are Ao = 1536 nm, 
r = 3.55 ns, L = 51 m, 7 = 3.26 W^km" 1 , (3 2 = -15.26 ps^m" 1 and Aft = 286.1 fs 
m _1 . The peak powers Pq corresponding to simulations (experiments) are (a) 132.5 W 
(125.6 W), (b) 127.2 W (121.4 W), and (c) 116.6 W (112.9 W). The numerical results have 
been convolved with a response function to take into account the resolution of the OSA. 
Amans5 7847F4. eps 
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